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We demonstrate theoretically that most of the observed transport 
properties of graphene sheets at zero magnetic field can be ex- 
plained by scattering from charged impurities. We find that, con- 
trary to common perception, these properties are not universal 
but depend on the concentration of charged impurities n^^p. For 
dirty samples (250 x lo'^cm^- < «i„p < 400 x lo'^cm^-), the value of 
the minimum conductivity at low carrier density is indeed 4e-/h in 
agreement with early experiments, with weak dependence on im- 
purity concentration. For cleaner samples, we predict that the min- 
imum conductivity depends strongly on «i„,p, increasing to ?,e~/h 
forni„,p ~ 20x lO'" cm ~. A clear strategy to improve graphene mo- 
bility is to eliminate charged impurities or use a substrate with a 
larger dielectric constant. 

graphene | electron transport | minimum conductivity 

The past two years have seen a prohferation of theoretical and 
experimental interest in graphene. The interest stems mainly 
from the striking differences between graphene and other more well 
known semiconductor-based two dimensional (2D) systems, that 
arise mostly from its unique band structure, obtained by consider- 
ing graphene to be a single sheet of carbon atoms arranged in a 
honey-comb lattice. Graphene is in fact a carbon nanotube rolled 
out into a single 2D sheet, and as was already known from studying 
carbon nanotubes, electrons moving in the periodic potential gener- 
ated by the carbon lattice form a band that displays striking properties 
such as having a chiral Dirac equation of motion with a mathemati- 
cal structure similar to Weyl neutrinos. This intriguing 'relativistic 
Dirac-Weyl spectrum of graphene has attracted substantial interest 
and attention. While this analogy of considering graphene as a solid- 
state realization of the "massless chiral Dime Fermion" model (de- 
veloped as a solution to Dirac's Lorentz invariant generalization of 
Schrodinger's equation) has some utility, it also has the potential to 
be misleading. In particular, as we argue here, searching to explain 
the experimental transport properties of graphene by focusing on the 
"Dirac point" (see formal definition below), obscures the real mech- 
anism of carrier transport. In this respect, and from our perspective, 
the physics of graphene has more in common with the Metal-Oxide- 
Semiconductor-Field-Effect-Transistors (MOSFETs) that form the 
backbone of our current day semiconductor industry, than with the 
physics of relativistic chiral Fermions. Studying graphene is there- 
fore as much about making useful MOSFETs from pencil smudges, 
as it is about studying quantum electrodynamics in a pencil mark. 

We observe that already within one year since the fabrication of 
the first gated 2D graphene samples (that enable a variable external 
gate voltage tuned carrier density), mobilities as high as 2.5m^/Vs 
have been reported, and these values are comparable to the best Si 
MOSFET samples at low temperature. In addition, graphene mobil- 
ity is relatively temperature independent, making room temperature 
2D graphene mobilities to be among the highest in FET type devices. 
It is therefore both of fundamental and technological interest to un- 
derstand the transport mechanism in graphene in reasonable qualita- 
tive and quantitative detail. Similar to MOSFETs, transport proper- 
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Fig. 1 . (A) Main panel shows cartoon of our model where charged impurities 
in the substrate at a distance d away from the graphene sheet create a spatially 
inhomogeneous screened Coulomb potential. At low carrier density the system 
breaks up into puddles of electrons and holes. The residual density is then cal- 
culated self-consistently giving remarkable agreement with experimental results. 
The inset (taken from 5 ) is used here just to illustrate the voltage fluctuations 
schematically, where it is understood that the presence of both electron and hole 
carriers implies that both positive and negative voltages are screened. (B) Pre- 
dicted conductivity traces (Eq. [U for different values of /limp- Curves are offset 
vertically by lOOe-//! for clarity and show from top to bottom, impurity concentra- 
tions (in units of lO'^cm"^): 20 (very clean), 40, 80, 160, 320 (very dirty). Most 
of the samples considered in Ref. 6 had mobilities in between the bottom two 
curves. Clearly seen are the gate voltage offset and the minimum conductivity 
plateau, which are larger for the dirtier samples. 

ties of graphene are determined by scattering from charged impurities 
which are invariably present. We report here the essential graphene 
transport theory focusing on charged impurity scattering. 

Prior to this work, the conventional wisdom in the graphene com- 
munity was that close to the Dirac point, carrier transport is ballistic, 
that the minimum conductivity is universal and that we lack a basic 
understanding of how at the Dirac point there could be carrier free 
transport over micron-sized distances. In this context, our work pro- 
vides a theoretically simple explanation for this graphene transport 
mystery: charged impurities in the substrate generate carrier density 
fluctuations that allow for non-universal diffusive transport, and that 
these density inhomogeneities render the Dirac point physics experi- 
mentally inaccessible, at least for current graphene samples (see note 
added at the end of this paper). 

We emphasize that although the importance of charged impurity 
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scattering in determining tlie linear-in-density high-density (i.e. far 
from the charge neutral Dirac point) graphene carrier transport has 
already been established (Tl[2l[3l|4l in the literature, the current work 
is the first to point out that the same charged impurities will have a 
qualitative effect at low carrier density close to the Dirac point by 
providing an inhomogeneous electron-hole puddle landscape where 
the conductivity will be approximately a constant over a finite range 
of external gate voltage, providing a simple and physically appeal- 
ing explanation for the observed graphene minimum conductivity 
plateau. The importance of our work therefore lies in its ability to 
explain the graphene transport data at low carrier density through 
a physically appealing charged impurity-induced mechanism which 
quantitatively explains both the existence (and the width) of the min- 
imum conductivity plateau as well as its magnitude. No such expla- 
nation existed in the literature before our work, and the minimum 
conductivity phenomenon in graphene was considered to be an out- 
standing experimental puzzle. 

In this context, ''bare graphene" is the empty honeycomb lat- 
tice, where allowing for electron hoping between adjacent sites gives 
the linear Dirac-Weyl spectrum. The nominally undoped or ungated 
situation is that of the completely filled valence band and a com- 
pletely empty conduction band which touch at the Dirac point mak- 
ing graphene a zero-gap semiconductor. The Dirac point is a singular 
point of measure zero that separates the conduction and the valence 
band in the linear graphene spectrum. This "intrinsic graphene" with 
the chemical potential (or Fermi energy) precisely at the Dirac point 
has no free carriers and is obviously an abstract model, since the 
slightest amount of doping or external potential will induce carriers 
in the system. Charged impurity disorder or spatial inhomogeneity 
will render this intrinsic graphene experimentally unrealizable. "Ex- 
trinsic graphene" is when one induces free carriers (either electrons 
in the conduction band with a positive potential, or holes in the va- 
lence band with a negative potential), by applying an external gate 
voltage Vg, or equivalently, by doping the system. Carrier bands are 
filled up to a certain Fermi Energy Ep, determined by the electrostatic 
potential configuration of the graphene environment, and depending 
on which is larger, could be dominated either by the external gate 
voltage or by charged impurities. All experimental graphene samples 
are extrinsic, since there are invariably some free carriers present in 
the system, and transport close to the Dirac point is dominated by 
two distinct effects of the charged impurities in the system: (i) the 
induced graphene carrier density is self-consistently determined by 
the screened charged impurity potential; and (ii) the conductivity is 
determined by charged impurity scattering. 

As one would expect with any newly discovered electronic ma- 
terial (7), the full gamut of experimental techniques has been used 
to explore graphene properties, including light scattering |8|, angle 
resolved photo-emission spectroscopy |9|, and surface probe mea- 
surements 1101 II II . The field continues to evolve with new ideas 
being explored including making suspended graphene fT2l| and elec- 
tromechanical resonators 1 13|, transfer-printing graphene onto plas- 
tic [14 1, using superconducting |15| and ferromagnetic leads |16|, 
having patterned top-gates [17. ,18. .191, exposing samples to molec- 
ular dopants 1201 and fabrication of graphene-nano-ribbons 1211 . 
These rapid experimental advances show that the study of graphene 
is still in its infancy, with many promises for the discovery of new 
physics and for application to technology. 

The focus of the present work is on the important graphene trans- 
port measurements l6l 1221 . These were the first experiments to be 
done, and are still not understood, in that there are widely disparate 
claims on the transport mechanism. Since these experiments form 



the basis for most of the future work on graphene as well as its 
prospective technological applications, a correct understanding of the 
basic transport physics is of fundamental importance. While there are 
many features observed in the transport experiments, two have been 
highlighted and particularly discussed in the literature, namely, the 
low-density "minimum conductivity" ctq, i.e. the value of the conduc- 
tivity at or near the Dirac point (where Ep ~ 0, and a naive picture 
would suggest that there are no charge carriers); and the high-density 
conductivity a(n) which is linear in the carrier density n, giving a 
constant mobility p = a/ne. As we discuss below, our analytic the- 
ory is the first to explain both features quantitatively, as well as make 
predictions for the width of the minimum conductivity plateau and 
the offset of the Dirac point from zero gate voltage. 

In the literature, most theoretical work has focused on the short- 
range scattering mechanism (also called "white noise" disorder) 
to understand graphene transport, mainly as a matter of technical 
convenience. Early work using a Kubo formalism in the ballis- 
tic limit 1231 1241 showed that the conductivity for massless Dirac 
Fermions is e-/(nh) for vanishing disorder, and that this universal 
value occurs only at the Dirac point and not in its vicinity. At fi- 
nite carrier density, the Kubo formalism with short-range scattering 
gives a conductivity that is constant with carrier density (Tj and not 
the linear in density behavior seen in experiments. Certain numer- 
ical |2| and analytical |25| methods that try to extrapolate between 
these two limits inevitably get a square-root dependence of conduc- 
tivity on density (not linear), and give orders-of-magnitude incorrect 
values for the mobility. More recently, short range scattering has 
been considered theoretically 1261 |27l 1281 with the finding that at 
zero temperature, localization effects should give ctq = 0. While 
all these works improve our abstract theoretical understanding of 
graphene, they are all in qualitative disagreement with existing ex- 
perimental data. We argue here that short-range scattering has little 
to do with the experiments of Refs. |7, 6, 22 1, and that although 
localization effects may very well be important in the zero temper- 
ature limit, the existing bulk graphene transport data at accessible 
temperatures (T > OAK) are in the Drude-Boltzmann diffusive trans- 
port regime. Equally important, the observed transport properties of 
doped graphene do not access the Dirac point physics, at least in the 
currently available samples |6, 22 1 which have fairly large concen- 
trations of charged impurity centers. While the observed conductiv- 
ity value (and not conductance) of 4e^/h in dirty samples brings to 
mind connections with universal conductance quantization phenom- 
ena such as one-dimensional point-contact conductance quantization 
or Zitterbewegung or quantized Hall resistance, we argue here that 
no such universal physics is at play in current bulk graphene trans- 
port experiments where conductivity, and not conductance is being 
discussed. 

Our goal here is to develop a quantitatively accurate analytic 
theory for the most important graphene transport problem, namely, 
the regular bulk dc conductivity studied in the diffusive Drude- 
Boltzmann limit, which is the regime of technological interest. We 
argue that charged impurity scattering is responsible for most of the 
observed bulk diffusive transport behavior in graphene. Charged im- 
purities could reside either inside the substrate or created near the 
graphene-substrate interface during the processing and handling of 
samples. The typical concentration of charged impurities in a Si02 
substrate is «i,np ~ 50 x lO'^cm"-, and are known to dominate the 
transport properties of other extensively studied 2D semiconductor 
systems 1291 . 

Earlier work p2"3l'4l demonstrated that while the mean free path 
for short-range scatterers f j ~ 1 / ^/n, where n is the carrier density. 
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for Coulomb (i.e. charged impurity) scatterers ~ V"- This implies 
that although graphene at low density is a clean or ballistic system for 
short-range scatterers, it is a dirty and diffusive system for long-range 
scatterers. At low carrier densities of n ~ ni„,p ~ 50 x lO'^cm"^, one 
can estimate that > 1000 nm and tc ^ 50 nm. Therefore, at the 
lowest densities close to the Dirac point, graphene transport proper- 
ties are completely dominated by Coulomb scattering. This simple 
argument also establishes that any short-range scattering, even if it is 
present in graphene, is irrelevant in the low carrier density limit near 
the Dirac point where the minimum conductivity plateau is observed. 

In the present work we provide a complete picture of transport 
in both high and low carrier density regimes using a self-consistent 
RPA-Boltzmann formalism where the impurity scattering by the 
charged carriers themselves is treated self-consistently in the Ran- 
dom Phase Approximation (RPA), and the dc conductivity is calcu- 
lated in the Boltzmann kinetic theory. We derive analytic expressions 
for (i) mobility /j, (ii) plateau width, (iii) minimum conductivity, and 
(iv) shift in gate voltage. Our results for graphene on a SiOa substrate 
can be summarized as 
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where Eq. [TO]below gives analytic expressions for n* and h. 

An important finding of our work is that for dirty samples with 
"imp ~ 3.5 X lO'^cm"-, (To is indeed close to 4e^//i as observed ex- 
perimentally, and not very sensitive to changes in disorder, while for 
cleaner samples with n^^^ ~ 2 x 10" cm"-, ctq ~ ie^/h, and is sen- 
sitive to the value of «i„,p, thus explaining the mystery of why more 
recent experiments show a larger magnitude and larger spread in the 
value of (To. For the typical densities used in the early graphene ex- 
periments, the minimum conductivity appears to saturate at a univer- 
sal value of 4e^/h, but we predict that there is nothing universal here, 
and for dirtier samples, the value of the minimum conductivity as a 
function of «j,„p would slowly decrease. 

The formalism developed here to obtain Eq. [T]can be divided into 
three steps. First, we develop an analytical solution for the Boltz- 
mann transport theory using the full RPA treatment of the charged 
impurity scattering. We find analytically that a « 20(e^//!)(«/nj,np) 
in agreement with earlier numerical calculations j2][3l|4). Second, 
we extend the methods of Refs. L5^ .30] to evaluate the screened volt- 
age fluctuations induced by charged impurities. Here we calculate 
the potential fluctuations using the full RPA screening, which al- 
though being complicated and cumbersome, is necessary to obtain 
quantitative agreement. Third, we develop a theory to calculate the 
residual carrier density n* self-consistently. We find that the ratio 
"*/"imp> that is directly related to the minimum conductivity through 
(To ~ 20(e^//j)(n*/«inip), is a monotonically decreasing function of 
nimp and the dependence gets weaker for larger impurity density. 

Before we provide details of our calculation, we first address 
the range of validity of our self-consistent RPA-Boltzmann theory. 
First, we consider only Coulomb scattering. As already discussed 
above (see also Ref. |4|), other scattering mechanisms are irrele- 
vant at low density and the experimentally observed linear depen- 
dence of conductivity with carrier density singles out charged im- 
purities as the dominant scattering mechanism. Only in the limit of 
very small charged impurity density, must one include short-range 
scattering into the formalism (such short-range scattering may arise 
from point-defects and dislocations in the lattice). We find tTo = 
(4e^//!)[nimp/(5«*) -l- where we estimate r] = IjkfCs < 1/10, 
suggesting that for very low impurity densities (two orders of magni- 



tude lower Himp than present day samples), (Tq will saturate at around 
20e^/h, before our charged impurity model gives way to short-range 
scattering. Second, our theoretical calculation is done at zero temper- 
ature. Theoretically one expects | 3 1 very weak temperature depen- 
dence for T <K Tp, and this is indeed consistent with experimental 
observations where T/Tp < 0.2. Third, since the sample sizes are 
several microns in length, and the mean free path is tens of nanome- 
ters, we are certainly in the diffusive as opposed to ballistic transport 
regime. Fourth, electron interactions are treated within the RPA ap- 
proximation scheme. RPA is an expansion in r,, and in graphene 
experiments on SIOt, x 0.8, so one would expect it to work better 
than for 2D semiconductors or metals, where although rj ~ 2 - 10, 
RPA provides an excellent approximation. Fifth, the scattering time 
is calculated using Boltzmann kinetic theory. Formally, Boltzmann 
theory (as described by Eq. |3]below) is valid for kpC » 1, but it is 
also the standard theory used to describe Coulomb scatterers in 2D 
systems I29II31I . For clean graphene samples, fcpf > 4 at low den- 
sity and kp{ ~ 100 at high density, making the Boltzmann theory 
valid. Sixth, the average carrier density is obtained from the potential 
fluctuations using a local density (also called Thomas-Fermi) approx- 
imation as is normally done in 2D systems l5l l32l . This is valid so 
long as n*/"inip ^ 0.01, which guarantees that the spacing between 
electrons is much less than the length of the conducting cluster [33 1. 
This condition holds both empirically and is established a posteriori 
to be valid throughout the Boltzmann transport regime. Similarly, 
this demonstrates that in current experiments the Dirac cone is al- 
ways filled with electrons or holes whose average density ranges from 
0.2 njinp to 0.5 Himp. Formally, one can consider the limit ni,np — > oo, 
and find that n' = (?iinip/32d^;r)''^, with no lower bound on the min- 
imum conductivity ctq ~ 1/ ^/ni^, but this "mean-field" approxima- 
tion breaks down when the inhomogeneities become so large that 
calculating the conductivity through the Boltzmann transport of the 
average density n' becomes meaningless. Seventh, we ignore An- 
derson localization. The observed experimental absence of Anderson 
localization is not unique to graphene. In 2D semiconductor sys- 
tems 1341 1351 [36l . a percolation metal-insulator transition is seen at 
low carrier densities instead of a localization transition, whereas in 
graphene the percolation transition is an electron-metal to hole-metal 
transition without any intervening insulating phase. Moreover, recent 
numerical work (371 [381 suggests that there may be not be a metal- 
insulator transition in graphene. Naturally, one could relax any of 
the above assumptions in a future calculation, but for the purpose 
of comparison with current graphene experiments, these approxima- 
tions are well justified and provide a consistent framework to under- 
stand graphene transport. We observe that in our RPA-Boltzmann 
theory, (Tq reaches the so-called universal Dirac point minimum con- 
ductivity value I23II24I I39I of e-/(nK) for unphysically large impurity 
densities of nin,p > lO'^^cm"^, where the diffusive transport approxi- 
mation no longer applies. 

The theoretical picture presented here is shown heuristically in 
Fig. [U Charged impurities either in the substrate or in the vicinity 
of graphene create a spatially inhomogeneous potential distribution 
in the graphene plane. At low carrier density, the spatially inhomo- 
geneous potential breaks the system up into puddles of electrons and 
holes. This theoretical prediction has now been verified in a re- 
cent surface probe experiment |T 1 1 using a scanning single-electron- 
transistor to directly measure the potential fluctuations in graphene, 
and finds quantitative agreement with earlier predictions 1 30] for the 
height and width of the electron and hole puddles. In addition, there 
is recent indirect experimental support f 18', T9\ for the electron-hole 
puddle picture proposed in Ref. L4J. Unlike usual 2D systems, both 
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electrons and holes screen the external potential. These potential 
fluctuations directly change the local chemical potential inducing a 
residual density which in turn changes the screening. Here we use a 
self-consistent procedure to determine the residual density n' , which 
manifests itself in experiments by a residual conductivity plateau that 
is shifted by an offset gate voltage V° = hja, whose width is n* and 
whose magnitude is (t(«*), where a{n) is the RPA-Boltzmann con- 
ductivity for carrier density n. Vg is called the Dirac gate voltage 
because it is the value of the external gate voltage where the Hall co- 
efficient changes sign indicating that carriers change from electrons 
to holes. In MOSFETs, the voltage corresponding to is often re- 
ferred to as the '^threshold voltage", since in these systems it marks 
the onset of conductivity at this critical value of carrier density, and 
cr = below this threshold. In graphene, V" separates the conducting 
electron and conducting hole transport regimes, whereas in MOS- 
FETs the threshold voltage separates conducting and insulating 2D 
channels. Here a ~ 7.2 x lO'^cm^^V"' is a geometry related factor 
(that can be measured directly l6l l22l using Hall measurements) and 
is used to convert the experimentally measured gate voltage to car- 
rier density n. For the rest of the paper we develop the theory only 
for carrier density n = k^/n, where it is understood that any compar- 
ison with the experimentally measured gate voltages is made using 
n = aVg. 

Based on estimates from surface probe measurements (7] 1101 , 
and consistent with earlier work both in graphene [4j .30J and in Si- 
MOSFETs 1291 1401 , we assume that the charged impurities lie in a 
plane at a distance d ~ I nm from the graphene sheet and calcu- 
late the voltage fluctuations taking into account of screening using 
the RPA approximation. The screened voltage fluctuation is a func- 
tion of d and the carrier density n; a larger carrier density more ef- 
fectively screens the charged impurities, while the potential fluctua- 
tions are larger for low carrier density. We include this effect self- 
consistently in our theory, where n is both determined by and de- 
termines the screened impurity potential. Our theoretical results do 
not depend in any qualitative manner on the precise choice of d, and 
one can develop relationships between the four experimental quanti- 
ties (i.e. mobility, plateau width, minimum conductivity and shift in 
gate voltage) that are independent of d. We emphasize that since a 
single parameter iii^p determines all four experimental quantities, we 
anticipate that our theory would be consistent with each of them to 
within a factor of 2. Comparison with representative samples from 
the Columbia, Manchester, and Maryland groups (see inset of Fig. [2] 
and note added at the end of this paper) shows agreement for the mo- 
bility and gate- voltage shift and agrees to within the expected factor 
of 2 with measurements of the plateau width and minimum conduc- 
tivity. 

We now proceed to calculate the Boltzmann transport conduc- 
tivity. For 2D graphene, the semi-classical diffusive conductivity is 
given by 
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where = g„ = 1 are the spin and valley degeneracy factors and 
the mean free path £ = vpT, with the scattering time t being given at 
r = Oby 
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where V(q) = Ine^'''' I (Kq) is the Fourier Transform of bare 
Coulomb potential at the transfer momentum q = |k - k'| = 
Ikf sin(6/2). While the exact RPA dielectric function is known 1411 . 
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Fig. 2. The main panel shows Eq.[To] The two lines starting at the origin show 

"/"imp for ;iimp = 10 X lO'^cm"^ (very clean) and 

'^imp — 350 X 10 cm 
(very dirty). The RPA and Thomas-Fermi results show 2r^Co{rs,a = 4d^/mi) 
as a function of carrier density n. The points of intersection represent the self- 
consistent solution, and one can read off n' from the x-axis, ;i*/;iimp from the left- 
axis and the "minimum conductivity" o"o from the right axis. The inset shows o"o 
as a function of charged impurity concentration showing that (i) it is non-universal, 
(ii) dirty-samples have o-q = 4e^//i over a wide range of impurity concentration 
and (ill) the cleanest samples have o-q ~ 8e-//!, but the value is sensitive to 
the concentration of charged impurities. Also shown is comparison with repre- 
sentative experimental results, where the square shows results from Columbia, 
diamond from Manchester, and circle from Maryland. These same three sam- 
ples (with conductivity shown over the full density range) were compared to a high 
density numerical Boltzmann theory in Ref. 



for our purposes we can use the following simple approximate ex- 
pression which allows for analytic calculations and provides results 
that are indistinguishable from the exact results (see Fig [3} 
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Fig. 3. Comparison of voltage fluctuation CoUs.a) using different screening ap- 
proximations. The Random-Phase-Approximation (RPA) is the main approxima- 
tion used in the present work. The Thomas-Fermi result was derived in Ref. [30] 
and the "complete screening" result valid for rjO » 1, was obtained by Ref. "5]. 
All three approximations agree in the large density limit, but disagree for small 
density. Shown in dashed lines are small density analytic asymptotes for the 
Thomas-Fermi and RPA and squares show the numerical evaluation of Eq. |8] 
using the exact dielectric function reported in Ref. RTl . 
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where q,, = 4kprs, and solving the integrals exactly we find 
e- n 1 
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where for graphene on a Si02 substrate, = e^/(hyK) « 0.8, 
G(2r,) X 1/10 and tr a 20{e^ /h)(n/n,„^p). Note that G[x] is positive 
and real for all x. From this we derive a simple analytical expression 
for linear tail mobility in graphene 

ii.50^, 

/^O "imp 

where /jq = lm^/\^i and no = lO'^cm"^. Note that Eq. |6]depends 
only on the charged impurity scattering concentration iij^p, indicat- 
ing that the only way to increase graphene mobility for fixed r, is to 
improve the sample quality. However, Eq. [5]also shows that mobility 
depends on the substrate dielectric constant k, and therefore chang- 
ing the underlying substrate from Si02 to one with a higher dielectric 
constant would reduce r„ and would be another way to increase sam- 
ple mobility. For example, changing the substrate to Hf02 (k^ ~ 25) 
from Si02 (/c, ~ 4), should enhance graphene mobility by a factor of 
5. 

We now calculate the statistics of the random voltage fluctuations 
significantly extending earlier numerical work [5 1 to incorporate ana- 
lytically the non-linear screening of electrons and holes in a zero-gap 
situation, to find 
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where the superscript on C^^^ indicates that we used the RPA Ap- 
proximation, and El (z) = f_ r^e^'dt is the exponential integral func- 
tion. The voltage fluctuation result Coir^ = 0.8, a) is shown in Fig. [3] 
comparing different approximation schemes used in the literature. 
The analytic result Eq. [8] is compared with a numerical evaluation 
of Eq. |7]using the exact dielectric function first reported in Ref. 1411 . 
as well as with the long-wavelength (also known as Thomas-Fermi) 
approximation l30l , where e(q) - l+qJqfoiiLllq, and the "complete 
screening^' approximation f2','5l that is valid only for r^a » 1, where 
e(q) = qjq. In the limit of large a (i.e. kpd » 1), both the Thomas- 
Fermi and the RPA result approach the complete screening limit of 
C^^(rs,a) = {2rsaY~, that was previously obtained in Ref. (5|. The 
three approximations disagree in the small a (or low density) limit 
where 



0): 



-1 



In 



2r, 



2r, + 1 



4 ln(7a) 
(2 + Tir.f 



[9] 



2r, + 1 

and 7 a 1.781 is Euler's constant. Notice that for RPA, a combi- 
nation of the small and large a = 4fcpd asymptotes span most of the 
density range. While these analytical asymptotes correctly describe 
the screened potential for most densities, it turns out, that the regime 
relevant to graphene experiments is the window where they do not 
work well, and the full functional form of Cq^^ shown in Eq. [8]needs 
to be used. 

As discussed earlier, to determine the self-consistent residual 
density, we equate the average chemical potential to the fluctuation 



in the screened charged impurity induced potential as E^ = 6V-, and 
find 



n 



— = 2rJCo^PA(r„fl = 4rfV^), 



n 



imp 

4n* ' 



[10] 



where the second expression for the impurity induced shift in volt- 
age is determined from V = OTii„,p7/(2fcp) |5] |30I . Combining these 
results gives Eq[T] 

Shown in the main panel of Fig. [2] are the results of the self- 
consistent procedure. The inset shows the value of the minimum 
conductivity and compares with the experimental results. These are 
the same three samples that were shown in Ref. |4J to compare the 
conductivity at high carrier density (far from the Dirac point) with a 
numerical Boltzmann theory, and here we show the low carrier den- 
sity comparison near the Dirac point. Through Eq. [S] we have the 
high-density measurements directly giving )ii„,p, and this is the only 
parameter used to determine the minimum conductivity ctq. Our re- 
sults show that contrary to common perception, the graphene min- 
imum conductivity is not universal, but that future cleaner samples 
will have higher values of cro. We emphasize that, in addition to 
explaining the value of cro and its dependence on the sample qual- 
ity, our theory also naturally accounts for the width of the minimum 
conductivity plateau in agreement with experiments. For example, 
Kimp = 350 X lO'^cm"^ gives n' = 70 x 10"'cm"- and plateau width 

AVg = lOV. 

We emphasize that the most important qualitative result of our 
theory is to introduce a realistic mechanism operational in all dis- 
ordered graphene samples (i.e. in the presence of random charged 
impurities) which produce a plateau-like approximate non-universal 
minimum graphene conductivity at low induced carrier density. We 
obviously can not rule out other possible "universal mechanisms" 
which will lead to a "universal" minimum intrinsic graphene con- 
ductivity at the Dirac point in the clean limit, a situation beyond the 
scope of our theory. But the fact that the currently existing exper- 
imental data from three different groups exhibit non-universal min- 
imum conductivity in approximate (within a factor of 2) agreement 
with our theory indicates that any intrinsic universal mechanism be- 
yond our model may not yet be playing any role. We note however, 
that there have been recent theoretical I39ll42l|43l and experimental 
1441 work dealing with ballistic transport in mesoscopic graphene 
which show universal behavior in the regime where £ > W > L, 
where { is the mean free path, and W and L are the sample width and 
length. Our theory does not apply in this zero disorder ballistic limit 
since our work is entirely built on the picture of diffusive transport 
through disorder induced electron-hole puddles. 

In summary, we believe we have qualitatively and semi- 
quantitatively solved one of the main transport puzzles in graphene, 
namely, why the experimentalists see a conductivity minimum 
plateau and the extent to which this minimum conductivity is or is 
not universal. The theory developed here should only be taken as 
the first step toward a full quantitative theory of graphene transport, 
particularly at the lowest carrier densities. Many questions still re- 
main open, although we believe that we have taken an important step 
in the right direction. In particular, the precise nature of transport 
at the charge neutral Dirac point can not be accessed by our self- 
consistent treatment which is valid only at finite doping away from 
charge neutrality. The percolation through electron and hole puddles 
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becomes precisely equivalent at the Dirac point, and a purely percola- 
tive theory, not our self-consistent RPA-Boltzmann theory, would be 
necessary to understand the transport. We do, however, mention that 
the observed smooth behavior of conductivity as a function of gate 
voltage through the charge neutrality point indicates a lack of any 
dramatic phenomena at the Dirac point, and given that our theory is 
a good description of transport away from the Dirac point, it is con- 
ceivable that it remains quantitatively valid at the charge neutral point 
also. Given the great deal of current interest in graphene and the 
potential for graphene-based electronics applications, our transport 
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